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ABSTRACT 


An algorithm was designed to load different sized rect- 
angular solids into a container. It allows the option of 
forming pallets of material before loading the container. 
The algorithm will permit loading of cargo that may or may 
not be used as load bearing support for other cargo. Cargo 
1S allowed to be rotated if desired to improve efficiency 
and both the pallets and the shipping container may contain 
"voids" or volumes in which cargo is not permitted. A test 
of the algorithm utilizing an actual cargo list showed two- 
dimension (area) efficiencies of 95% and three-dimension 


(volume) efficiencies of 89%. 
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I. INTRODUCTION 


A. BACKGROUND 

Computerized analyses and computer assisted algorithms 
have been utilized extensively in most areas of transporta- 
tion systems. The military, in particular, has relied 
heavily on loading simulations and computer assisted algo- 
rithms to predict the assets required to meet a given trans- 
portation demand [2]. However, no reference can be found 
which indicates that these computer techniques have been 
accurate, flexible, or descriptive enough to act as an actual 
blueprint for loading multicommodity cargo into the trans- 
portation container, be it a sea van, truck or airplane. 

The actual loading is apparently still performed, for the 
most part, by personnel without the help of computers. 

Biman attempt to partially fill this void, an heuristic 
algorithm was developed which should be efficient and precise 
enough to use as an actual blueprint for loading one, two, 
or three dimension cargo. This algorithm exceeded the un- 
assisted performance of loading crews for sample data; is 
adaptable to any shaped container; permits container 'voids" 
or volumes where cargo can not be loaded (i.e., refrigera- 
tion vests, reserved space, etc.); recognizes that some 
cargo may be rotated for increased efficiency while other 
cargo can not be rotated; permits only weight bearing cargo 


to be used as a base upon which to stack other cargo; and, 





allows the optional requirement that smaller boxes must be 
loaded on a standard pallet prior to being placed into the 
container (the formal term of placing cargo into the con- 
tainer is 'stuffing' as opposed to 'palletizing' the cargo 
ООШОТ το stuffing). In addition, the algorithm is capable 

of solving a large problem within several computer (CPU) 
minutes and requires relatively little main core memory. 

The algorithm described above is hereafter designated as 

the container stuffing algorithm. The need for such an algo- 


rithm is discussed in the next section. 


B. THE NEED FOR AN ALGORITHM 

The advent of mechanized warehouses, sharply increasing 
transportation costs, and increased availability of computers 
greatly increases the potential return on investment that is 
expected to be realized from the implementation of an algo- 
rithm as described above. 

Mechanized warehouses permit extremely rapid, efficient, 
and flexible issuance of material from the warehouse. 
Material in a mechanized warehouse is received, stored, and 
issued with very little manual intervention by the warehouse- 
man. This is accomplished by the use of real time data bases, 
one hundred percent visibility within the receipt-issuance 
cycle, and complete knowledge of the item characteristics 
(weights, dimensions, etc.) of the material being stored. 
However, some of the efficiency gained by the mechanized 
warehouse is lost once the material is dispatched from the 


mechanized warehouse. For example, the material must be 





staged in the shipping section prior to the actual loading 
of the material into a shipment container. This is neces- 
sary because of the current inability to accurately predict 
necessary transportation assets and because of the need by 
the loading personnel to physically view and study the 
physical characteristics of the cargo prior to commencement 
of the loading process. By eliminating the need for stag- 
ing material prior to shipment, savings could be realized 
in manpower necessary to actually load the material, in 
staging cost, and in costs associated with positioning the 
container prior to commencement of the loading. Additional- 
ly, there is a cost associated with delaying release of the 
container pending completion of documentation which must be 
prepared after the cargo is loaded but prior to releasing 
prercontalner. 

Considerable savings could also be realized if greater 
efficiency of cargo volume to container volume were possible. 
As one specific example, if Naval Supply Center, Oakland, 
California, could increase its efficiency for shipments to 
Japan and Philippines (705,400 cubic feet or 17,635 measure- 
ment tons per year) from the current rate of 80% to, say, 
87%, a yearly savings of transportation costs would be 
approximately $266,000 [12]. 

The algorithm developed to satisfy these requirements 
will be discussed in detail after the next section which 
briefly discusses the generic operations research problem 


which is becoming known as the loading problem [3]. 
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C. THE LOADING PROBLEM 

In order to understand the stuffing problem it is first 
necessary to review its superset, the loading problem. The 
generalized loading problem is one in which items, leI of 
magnitude q; and value σι, are placed in containers, с 
capacity E and cost of di. The sets I and C contain, re- 
spectively, all items to be loaded and all containers used 
in the loading. 


The problem may indicate: 


a. Br Lg, and all items are loaded, or 


јЕС var 
b. BGs E „а, and all items need not be loaded; or 
jeC Te T 
Le. < Σα... 
Πες ig I 


The objective may be to 
a. Minimize Σ (αι ον) where S is the set of all items 
TES 
not loaded (Sel); 


рат ое ὃ (c e d.) 
j G 


A matrix of the objective functions and problem state- 
ment is shown below and summarizes the various possible 
problems. An asterisk indicates that the assumption of 
additivity has been made. Under this assumption, whenever 
q, < = all i items may be loaded in the j containers; 
Lem е that is, the assumption of additivity allows the 


quantities which are being loaded to be added together with- 


out geometrical considerations of individual containers. 
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This assumption is easily made when measurements are in 


terms of money, weights, liquid volume, or when max(q,) << 
lg I 

ПОП с) апа prior palletization is not used. 

JEC 


Problem Statement 


Objective 





Figure 1. Loading Problem Subsets 


Problem ο is the classical multidimensional knapsack 
loading problem which has been extensively analyzed [8]. 
Problem 3. has been solved for the case where c 7 Cj TOT 
ντ .πη j by Eilon and Christofides [3]. The Problem 3 
solution for the case where items I are rectangular solids 
was developed by Gilmore and Gomory who used very large 
scale integer programming techniques [7]. The Problem 3 
solution where dj = dk for all i, k and items I were rect- 
angular solids with one set of common dimensions was given 
by Seam and Sivazlian [10]. Problem 4 solution where items 
I were rectangular solids was presented by DeSha [2]. The 
Problem 4 solution for the case where items I were parallele- 
pipeds and C was a single container was given by Galata and 
Stoyan [4]. This paper is concerned with a specific subset 
of Problem 4 called herein the container stuffing problem 


which is a representation of the generalized method typically 


T. 


used to ship cargo. Although the problem was formulated in 
terms of a shipping problem, it may be easily expanded to 


solve related problems such as those presented by Brown [1]. 


D. THE CONTAINER STUFFING PROBLEM 

This problem is one in which n boxes of size RE are to 
be loaded onto pallets of capacity py > which are, in turn, 
loaded (stuffed) into containers, Ci, ON capacity Cy: In 


P problem BOX. < p, < €, for all 1, к апа 7 вох. < Xp, < 

j — k — `i . j 
Lc. with an objective of minimizing the number of containers re- 
L 


quired to load a given series of . j=(1,n). Because of 


the relative closeness (in size) of Ὁ tO Py and Py tO οι, 
geometric considerations are extremely important in obtaining 
a feasible solution to the minimization problem; and, thus an 
elementary but important constraint must be addressed: none 
of the boxes (pallets, containers) may overlap into the 

space occupied by another box (pallet, container). 

A successful practicable solution to this problem must 
also consider these following points. A rectangular solid 
box n. be loaded into a container six different ways de- 
pending on the relative positioning of the moving coordinate 
system (x', y', z') associated with the box and the fixed 
coordinate system (x, y, Z) associated with the container. 

It is assumed the container has its "origin" located at 
o 00) with length, width, and height in the x, 
y, Z direction. Figure 2 shows the six possible orienta- 


tions of a box in a container. These degrees of freedom 


Τη 











Figure 2. Possible Orientation of a Box within a 
Pallet;Container. 
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permit n!6™ possible sequences for loading n. boxes. For 
instance, six boxes may be loaded into a container in more 
than 33 million different sequences. 

Because of this vast number of sequences, a practicable 
solution must be descriptive as well as prescriptive. It 
must define not only the sequence of the loading but also 
the relative position of each box in the Les and the 
relative position of each pallet in the container. 

Fortunately, certain real world constraints reduce the 
number of feasible solutions. The nature of some of the 
boxes requires that the box be loaded "this side up", that 
is, the box has a predetermined orientation (this reduces 
the problem to only two possible orientations). Also, the 
boxes have different load bearing capabilities resulting in 
the larger, heavier boxes being placed near the bottom of 
the stack. 

The loaded containers must meet specified maximum 
weights and distributed weight parameters and hence may not 
end up being completely filled. Also, some containers re- 
quire that voids be reserved to permit air circulation 
around vents, or to provide space for future cargo to be 
loaded elsewhere, etc. Finally, it is often the practice 
of loading personnel to utilize one of the larger boxes as 
the "pallet base' for equal size (length to length and width 


to width) and smaller boxes. 
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E. RELIANCE ON HEURISTICS 

No reference could be found that presents an exact solu- 
tion for the container stuffing problem. The complexity of 
the problem and need for rapid solutions exceed the capa- 
bilities of even the most sophisticated mathematical solu- 
tions currently available. Therefore, the solution techniques 
presented in this paper rely heavily on heuristics to approxi- 


mate the exact solution to the stuffing problem. 
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II. OBJECTIVES AND SCOPE 


The objective of this study was to develop a flexible 
m rj nm capable of Solving the above defined container 
stuffing problem. The algorithm which was developed pre- 
supposed that a computer would be required for the calcula- 
tions and that all necessary data on the items to be loaded 
were available. These necessary items include dimensions 
of the box and its load bearing capabilities. 

For tractability, each input box and the container was 
assumed to be a rectangular solid. This assumption as it 
related to the input boxes could be relaxed in actual prac- 
tice by defining a rectangular solid which superscribes the 
object to be loaded and by designating this rectangular 
solid as a non-load bearing box. Likewise, the assumption 
of rectangularity of the shipping container can be relaxed 
by defining a rectangular solid which is superscribed by 
the actual container and by defining voids within this rect- 
aneamear Selaa,. The accuracy of this approximation is simply 
a function of the scaling of the dimensions used in the 
alkoorithm. 

Weight distribution of items within the container was 
not explicitly addressed by the algorithm. Neither was 
center of gravity restrictions. However, these restrictions 
could be easily included by modification of peripheral logic 


Τη υπ α от си. 
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"ΠΗ STUFFING ALGORITHM 


A. GENERAL DESCRIPTION 

The stuffing algorithm (the FORTRAN program is contained 
in Appendix A) was designed to provide as much flexibility 
as possible, to be descriptive (i.e., describe the loading 
procedure in terms of relative positioning of each box)as 
well as prescriptive, and to run as fast as possible on a 
computer. As will be shown below, the amount of time re- 
quired to execute the algorithm is a function of the options 
used as well as the level of optimization desired. By vary- 
ing the input parameters the following options may be 
exercised: 

(1) Prior to stuffing the container, the input boxes 
may be first loaded onto a "Standard" pallet whose dimen- 
Sions are specified by input parameters. 

(2) Boxes larger than a certain size can be used as a 
base upon which to stack other boxes. 

(3) Before loading a box the algorithm may or may not 
require that the box be supported at each of the box's lower 
Tour Corners. If support is not required, boxes may “over- 
hang" or even be suspended with no support as may be desired 
when, say, designing an electronic component comprised of 
many subcomponents which may be held in place by wiring. 

(4) Boxes may be individually specified as non-load 
hearing boxes which denies their use as support for other 


boxes.  Non-load bearing boxes may still be overstacked 
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provided the overstacked box receives support from load bear- 
ims boxes. 

(5) The container is initially defined to be a rect- 
ντα Solid with rectangular solids (i.e., voids) cut from 
Boe Original rectangular solid. Thus, any reasonable geo- 
metric shape may be approximated. 

(6) Voids may be described either in pallets or in the 
container or both, provided the voids can be constructed of 
a series of rectangular solids. The void need not be con- 
tiguous to a boundary of the pallet or the container. 

(7) Voids placed in pallets and/or containers may either 
be load bearing or non-load bearing. 

(8) Boxes may be rotated from zero to five times in 
order to improve local optimization of the loading process. 

(9) Boxes may be svecified whereby no other box is 
allowed contiguous to one of the specified box's five remain- 
ing sides. 

СКОЈ е level of optimization, and, therefore, the 
amount of computer time is specified by input parameters. 

eum The algorithm may be used to load either three, 
two, or one dimensional objects (rectangular solids, rect- 


ειπες, ον 11465). 


B. “MEASURE OF EPFECTIVENESS 

The measure of effectiveness (MOE) used in the optimiza- 
tion sequence was defined as the total volume of input boxes 
divided by the volume of the containers into which the boxes 


were stuffed. The volume of the containers was defined as: 


E 





(N-1)(VOL) + (MAXW + CONW)(VOL) , 


where: N = Number of containers utilized 
VOL = Volume of one container 
MAXW = Maximum width utilized in the last container 


CONW = Container width 


This measure of effectiveness was devised in order: to 
penalize for any wasted volume on pallets as well as wasted 
volume in the container itself and to allow differentiation 
between various loading sequences for the case in which all 
boxes were stuffed into one container. 

As an example, the measure of effectiveness for the 
example problem solved in Chapter IV in Figure 3, is 0.733 
and is computed as follows. 

Total box volume is 228,096, container dimensions and 
volume are 60x60x96 and 345,600, respectively, and all boxes 
were stuffed into one container. 


The formula shown above then gives: 


MOE 7» (228,096):((1-1)(345,600)-(54:60)(345,600))20.733. 


C. PREVIEW OF THE ALGORITHM 

Before discussing the algorithm in detail, it is first 
necessary to describe the general approach to the stuffing 
algorithm and to set forth a few basic definitions. 

The algorithm loads one "pallet" (as defined below) at 
a time by inspecting each of the n boxes in the ordered 


input stream of boxes. Box number one is inspected first, 
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box number two is inspected second, and so forth until all 
n boxes have been inspected. If the box under inspection 
can be loaded without violation of one of the conditions 
(constraints) discussed in Section IIIA above, the box is 
loaded. If the current box does violate one of the con- 
straints, it is passed over and the next box in the sequence 
is inspected. If none of the boxes waiting to be loaded 
can be loaded, a new pallet is begun. Thus, the algorithm 
maintains "feasibility" while searching for one particular 
solution (1.e., a local minimum) to the stuffing problem. 
The loading procedure requires a decision to be made as 
to where an additional box may be placed. These locations 
(defined as "possible" origins) are limited by the algorithm 
to be the origin of the pallet or one of three positions 
relative to each of the boxes and voids previously loaded 
on the pallet as measured from the origin of the pallet. 
These positions are defined as the "x-position", "y-posi- 
tion", and '"z-position'" and are shown in Figure 3. The 
dimensions of the box being added are denoted by BOXL, BOXW, 
and BOXH corresponding to its length, width and height, 
respectively. A pallet which contains j boxes and voids 
will have (3j*1) possible origins. These positions were 
selected as possible origins because they limit the possible 
positions of the next box to a finite, manageable number of 
locations and a fast verification of feasibility. Finally, 
of all the "possible" origins, a subset of "permissible" 


Origins is defined. This subset of "permissible" origins 
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is determined by deleting from all possible origins those 
Grieins which have already been utilized by loaded boxes; 
by deleting the z-positions of all boxes defined to be non- 
load bearing boxes; by deleting positions at which none of 
the boxes still in the input stream can possibly fit; and 
by deleting positions at which it is desired to have no 
boxes contiguous to a loaded box's face (for example, if it 
be desired to have no box to the right (y-direction) of a 
given box, the y-position associated with the given box 
would not be included in the set of permissible origins. 
The order of inspection of those permissible origins is 
(1) tbe x-position ordered from the first loaded box to the 
last loaded box; (2) the y-position similarly ordered; (3) 
lastly the z-position likewise ordered from the first to 
the last loaded box. This order of inspection tends to fill 
the pallet in layers, always starting from the pallet's 
origin and progressing away from the origin in all directions. 
In order to determine if a box may be loaded at a given 
permissible origin, it is necessary to maintain a record of 
all previously loaded boxes and their relative positions in 
the pallet. This is ERES by maintaining a record 
of the preyiously loaded boxes' origins (defined as the 
"current" origins) and the boxes' dimensions in the x, y, 
and z directions. Thus, the first part of the feasibility 
question is answered by considering the box at a particular 
permissible origin and determining if the box is wholly con- 
tained within the space of the pallet and if the box does 


not intersect any previously loaded box or void. 
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The second part of the feasibility question must be ad- 
dressed whenever the problem input parameters require that 
each box must be supported. A box is considered to be sup- 
ported whenever all four of its lower corners rest upon a 
load-bearing box or void. This requirement, when exercised, 
does not allow any "overhang" of the box. 

A local minimum is obtained by attempting to move each 
box, aS it is loaded, toward the origin of the pallet. This 
is accomplished by determining if the box may be moved along 
one of the three directions, x, y or z, toward the pallet's 
origin. Movement is permitted only if the box does not in- 
tersect any previously loaded box or defined void. The box 
is moved in one direction at a time and movement is continued 
in an iterative fashion until no further movement toward the 
origin is possible. 

Finally, the term "pallet" is formally defined as a volume 
in which boxes are loaded. .AÀ "pallet" has dimensions of 
length, width, and maximum stacking height. A "pallet" does 
not, itself, occupy space. The dimensions of the pallet are 
determined by the input parameters. In the algorithm, two 
pascsxorepallets are used, a standard“ pallet and a "minimum" 
pallet. The "minimum" pallet actually defines the smallest 
sized box which will be allowed to serve as a pallet base. 

The algorithm operates by selecting the next box in the input 
Stream of boxes which is outsized to the minimum pallet. The 
dimensions of this box are then used to define the pallet base 


upon which to stack subsequent boxes. If no box is outsized 
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to the minimum pallet, the standard pallet is then used as 
the pallet upon which to stack boxes. A box is outsized if 
either the length of the box is larger than the length of the 
minimum or the width of the box is larger than the width of 
Giemminimum pallet. Thus, if prior palletization is not 
desired, the minimum pallet and the standard pallet are de- 
fined to be the same size as the container. Conversely, if 
it be desired to use the boxes themselves entirely to 
palletize the remaining boxes, the minimum pallet dimensions 
are defined to be zero. 

In order to clearly present the algorithm, the next sec- 
tion will briefly describe the notation used in the stuffing 
algorithm. Following this section, the algorithm will be 
stated. After the statement of the algorithm, a brief, 


Sample problem will be solved for illustration purposes. 


D. NOTATION AND DESCRIPTION OF MATRICES 

ΠΕ συ τπςσ [Orth the exact Stuffing algorithm, nota- 
tion will be briefly covered and the matrices used in the 
algorithm will be defined. 

There are n boxes to be loaded and their characteristics 
are contained in а matrix C. The original identification of 
an individual box is denoted as N. The subscript which refers 
to the order in the input sequence of box N is 0. Matrix C 
1s formally defined as the (px5) matrix of input boxes as 


follows: 
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N r4 BOXL, BOXW) BOXH; 


N r BOXL BOXW BOXH 
p p p 


where: N Boxzirdentitrearıon ot BOX, » SD 


Ө 
rg 7 number of boxes with dimensions of 
( BOXL, X BOXW , x BOXH,). 
Since there are n input boxes, L r, = n. Matrix C is used 


θ 
θ 
in the algorithm as a device to maintain a record of boxes 
yet to be loaded. 


The subscript which refers to the order in which box N 


is loaded onto the pallet is j. Thus No represents the 


J 


completed notation for box ordering. By N. = 25 is meant 
j 3 


that a box with identification number two was the fifth box 
in the ordered input stream and was the third box to be 


loaded. The voids, V, which may be defined on the pallets 


and the container are subscripted with an i. Thus Vi - V Ou S 


represent void one on the pallet currently being loaded 


whose dimensions are length, width, and height of VL, = VL, 
Lu VW. and УН; = VH... 


In order to permit a very rapid determination of whether 


respectively. 


a given box will fit at a given origin, a digital model of 
the pallet is established and updated with each box which is 
loaded onto the pallet. This digital model allows the 
determination of fit to be made through a series of very 


fast logic checks as described in the next section. This 
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model is defined as Array A which is an ((m+n)x7) matrix as 


follows: 
- SE ONE ^ (xq WE) Cy, VW, ) (z4 * VH, ) 
2 % Ya 2 (X5tVLp ) (YotVWo) (z9* VH, ) 
A = m A (x ЖУ) wen (an Hn 
N x y Z (x +BOXL ) (Yg +BOXW )(Z4 *BOXH. 
8+5 Omi mej mj mj mej mtj mej mej mj 
where: X;, Yi, 2; are the coordinates of the iu void 


whose length, width and height are VL,, VW,, and ҮН;, respectively, 
"ENS m. 
Nac is the identification number of BOX, whose 
length, width and height are BOXL, , BOXW, , BOXH, , eo n, 
a) «s n, and whose origin is located at coordinates 

"IAS 

Thus, matrix A column one identifies the box (or void), 
columns two through four identify the box's (or void's) 
location nearest the pallet origin, and columns five through 
seven in conjunction with columns two through four describe 
the volume occupied by the box (or void). 

To facilitate the selection of the next origin at which 
Peteco will attempt to load the current box, a logi- 


cal array of possible and permissible origins is established. 


By an extremely rapid scan of this model, defined as matrix B, 
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Memes tr orlelo 1S quickly determined. Matrix B is an 


((m#n)x3) matrix as follows: 


XORG, YORG, ZORG- 
XORG YORG ZORG 
m m m 
B = 
XORG, YORG, ZORG, 
. m+ j ° m+ j : me 
where: XORG, , YORG, , ZORG, , СОН л аге logical l variables 


which when true indicate that СЕ a VL, , yi); a * uer Zi); 


and (σι, y E + VHi)), respectively are permissible origins 


1” 1 


next to voids and XORG YORG ZORG SS n, are 


Bue SG Оше? 

ШОШ Са variables ΕΝ ΘΗΝ ЕСА indicate that 
UU Ca Opa op) Sos 76 BOW 2, 29 2; and 
(x. Š 7ρ > (uy + BOXH. )), a ету are elle 
is en to loaded ΑΝ For ease of notation, once a 
permissible origin has been selected for inspection, it will 
be designated merely as (ORX, ORY, ORZ). 

Each row of matrix B corresponds to a row in matrix A. 
The three elements in each row of matrix B correspond to the 
x-direction, y-direction and z-direction possible origins 
associated with each void and box on the current pallet (and 
described by matrix A). Of all the possible origins, the 
permissible origins are then defined by setting to a true 


value each element in matrix B which corresponds to the 


possible origin which is also a permissible origin. It is 
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precisely these permissible origins where attempts will be 


made to load additional boxes. 


E. EXACT STATEMENT OF THE STUFFING ALGORITHM 

To stuff n boxes (which are described by matrix C) into 
containers the following steps are used. 

1. Select the first ordered box which is outsized to 
the "minimum" pallet dimensions. Define the pallet base as 
this outsized box. If no outsized box is found, define the 
pallet as the standard pallet. 

2. Establish matrix A as the digital model of the pallet. 

3. Establish matrix B as the logical model of possible 
and permissible origins. 

4. Load the outsized box if one were found. Otherwise, 
load the first box from matrix C. Augment matrices A and B 
with an additional row to represent this box. Adjust matrix 
B as necessary to remove, if necessary, an origin from the 
set of permissible origins. 

S. Select the next box in matrix C. If no more boxes 
are left, go to step 12. 

6. Select the next permissible origin. The next per- 
missible origin corresponds to the next true element in 
column 1 of matrix B (x-position), followed by the next true 
element in the second column of matrix B (y-position), and 
finally followed by the next true element in the third column 
ορ ο ο. ο Zewosition), Call the selected origin (ORX, 


ORY, ORZ). If all origins have been tried and the box may 
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be rotated, go to step 11. If all positions have Dbeen tried 
and the box may not be rotated go to step 5. 

БО Determines if the box will fit at this origin. If 
it will not fit go to step 6 (The method of this determina- 
tion will be described below). 

7B. Improve the density of packing, if possible (This 
procedure will also be described below). 

8. If boxes must be supported (as discussed above), 
determine if the box is supported (as discussed below). If 
the box must be supported but is not supported, go to step 6. 

9. Load the box and augment matrices A and B with an 
additional row. Adjust matrix B to preclude any origin that 
may not be used. 

ШОО Go to step 5. 

11. Turn the box and go to step 6. 

12. If all boxes in matrix C are not yet loaded, go to 
step 5. If all boxes are loaded and the container has been 
stuffed, terminate the algorithm. If all boxes are loaded 
onto pallets but the pallets have not been stuffed into 
containers, move the pallets into array C, define the minimum 
and standard pallets as the container dimensions, and repeat 
the algorithm by returning to step l. 

Step YA is determined as follows. Inspect each pre- 
viously loaded BOX, , eo апат ејуесаи the se 


J 
of the current box BOX. because it would intersect BOX, ie 
K j 
A(j,2) « (ORX * BOXL, ) and A(j,5) » ORX and 
k 
A(3,3) < (ORY + BOXW, ) and A(j,6) > ORY and 
k 


A(3,4) < (ORZ + BOXH, ) and A(j,7) > ORZ 
k 
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where: A(j,k) is an element in matrix A and (ORX, ORY, ORZ) 
is the origin being inspected. 

If voids (V) were introduced, a similar inspection would 
be necessary of each V,, 0 < і < т whereby VL,, VW,, and VH; 


would be substituted for BOXL BOXW and BOXH 


Om UE DE 
g 3 J 
respectively. 
Step 7B is determined as follows. Inspect, one at a time, 
each possible direction of improvement (x, y, Z). Each 


direction of improvement is found by inspecting the origin 
(ORX, ORY, ORZ) under question and each row of matrix A. To 
determine if improvement be possible in the x direction, the 


following logic check is made on each row, k, of matrix A: 


A(j,2) » (ORX+BOXL, ) 
k 


ACj,4) > (ORZ+BOXH, ) or A(j,6) < ORY or 
k 


or AC],3) > CORY +BOXW , Jae Te 
k 


A(j,7) < ORZ. 


A true condition indicates that improvement is not possible 
at this row in matrix A. A false condition indicates that 
improvement is possible. The magnitude of improvement is: 

ORX - A(j,9) if improvement is possible or 99,999 if improve- 
ment is not possible. Now denote as slack, the magnitude of 
improvement found by inspecting row k of matrix A. The improve- 


ment found oyer all rows of matrix A is then: 


min(slack -Slack a ους. 


je +n 


To determine improvement in the y direction the logic 


check is: 


oul 





A(j,3) > (ORY+BOXW, ) or A(j,2) > (ORX+BOXL or 
k 


) 
К 


A(j,4) > (ORZ+BOXH or A(j,5) < ORX or 


) 
% 
A(j,7) < ORZ. 
The magnitude of improvement when the logic check is 
false is ORY - A(j,6). 
To determine improvement in the z direction the logic 
check is: 
A(j,4) > (ORZ+BOXH, IO AZ) > ( ORX*BOXL, or 
K k 
OO > (ORY*BOXW, MI AD) S ORX or 
k 


Eug 6) < ORY. 


The magnitude of improvement when the logic check is 
false is ORZ - A(j,7). 

With each improvement, the origin (ORX, ORY, ORZ) of 
the box being loaded is adjusted. The search for improve- 
ment is continued until no improvement is found in any of 


the three dimensions. 


Step 8 is determined as follows. Accept BOX, as sup- 
k 
ported if some BOX, "xL motor which Zo S BORE, = 
J J J 
Zo mies j < kas mtn sand BOX, is not declared a non-load 
k J 

bearing box, satisfies the following condition: 

CO < XX, апа =A 5) 2 XX, and 

ACJ,3) < yy, апа А(ј,6) > уу, 
for u = 1,2,3,4 where XX and yy, are defined as follows: 
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ll 
p— 


Ome Ὁ 


x 


When: u 


u gi 
J 
u = 2 or 4 Xx, = (Χα. i BOXL, ) 
j k 
и = lor 2 уун = Ye. 
J 
u = 3 or 4 uc ET * — 


Note the four (u) conditions are independently considered and 
one or.mare BOX, must b2 required to satisfy all these conditions. 
If load bearing voids were present, a Similar inspection of Vi, 


O < i < m, would be necessary. 


F. SAMPLE TWO-DIMENSIONAL PROBLEM 

In order to illustrate the algorithm, the following example 
is presented. It is a two dimensional problem since all features 
of the algorithm can he covered in two dimensions and two 
dimensions are more easily demonstrated. The three-dimensional 
stuffing algorithm is converted to a two-dimensional one by 
defining the height of each input box to be the height of the 
container, so that no stacking in the z direction occurs. 

In this example problem, 4 boxes (of base sizes 30x30, 
34x24, 20x20, and 26x10 and identification numbers 1 through 
4 respectively) are to be loaded without standard pallets and 
vuot prior palletization into a container of 60(L) x 60(W) 
x 96(H) with one yoid (10x10) which is located in the lower 
meat Geruer= ot the container.  For"illustration," "this "void 
may not have a contiguous box to its upper side. Each box 
may be rotated only once (i.e., all the boxes are marked "this 


side up"). The input stream of boxes has been ordered by area 


ο 








(length times width). Since only two-dimensional loading is 
considered, the boxes' heights are defined to be the height 
of the container (i.e., 96). The minimum pallet and the 
Standard pallet widths, lengths, and heights are 60, 60, and 


96, respectively. Matrix C has elements as follows: 


1 1 30 30 96 


2 1 24 34 96 


С = 
3 Jl 20 20 96 
4 l 26 10 96 

Step 1. Select the first box outsized to the minimum 


pallet. There were none, so define the current pallet dimen- 
sions to be 60x60x96. 

Step 2. Establish the A matrix which now contains all 
the voids (in this example there is only one). The A matrix, 


oor point is: 


A = (1 50 O Q 60 10 96 ) 


Step 3. Establish the B matrix of origins. Since the 
void has been specified so that no box may touch its upper 
tace its y-direction origin is not a permissible origin. 
The void also does not have a permissible x-direction origin 
nor z-direction permissible origin because no boxes yet to 
be loaded can fit at either of these origins inasmuch as 
the distance from these origins to a boundary of the pallet 
is zero. In fact, since all boxes have the same height as 
the pallet, B matrix will never show a permissible origin in 


the z-direction. Thus, B matrix is now: 
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БЕ КОЕ F) 


Ep Load the next box in matrix C. This is box one 
which has dimensions of 30x30x96. Augment the A matrix to 


include this box as follows: 


l O O O 30 30 26 


Augment B matrix and show permissible origins by setting the 


applicable element to true. Thus, 


Step 5. Select the second box in the C matrix (24x34x96). 

Step 6. Select the first permissible origin. That is, 
ο πο τις B Column by column always starting at the top of 
each column and working down. In this case, element B(2,1) 
is the next permissible origin. This element translates 
Eu origin of (A(2,5), A(2,3), A(2,4)) or (30, O, 0). 
Denote this as the current (ORX, ORY, ORZ). 

Step 7. Determine if the box will fit at this origin. 
This is accomplished by the following logic checks of matrix 
LE ocrucwscondition for any row of matrix A indicates that 
the box will not fit. Thus for, say, row one the check pro- 


ceeds as follows: 
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30 


20 


10 


BOX 2 


(24x34) 


BOBS 1 


(80x30) 


10 20 


BOX 4 
(26х10) 
ВОХ 3 
(20x20) 
| 
VOID | 
(10x10) 
30 40 90 60 


Figure 4. Example Problem. 


36 





AMS (ORXFBOXL ) and A(1,5) > ORX and 
2 
9 ) and A(1,6) > ORY and 


2 
2 


2 
A(1,3) < (ORY+BOXW 


A(1,4) < (ORZ+BOXH, ) and A(1,7) > ORZ 


This equates to: 


ОГ Оза апа 60 > 30 and 
O < ( 0+24) and 60 > O and 


ο ο ο ο) and 96 > O 


which is obviously true. Therefore, the box will not fit 
at this origin because it intersects the void. 
Step 6. Select the next permissible origin which is related to 
e (A(2,2), A(2,6), A(2,4)) or (0, 30, 0). 
Step 7A. Determine if box will fit at this origin. 
The following logic checks are made for rows one and two, 


respectively: 


50 < 34 and 60 > and O < 54 and 10 > 30 and 


O < 96 and 96 > and 


О < 34 and 30 > 


O 


and O < 52 and 30 < 30 and 


OT об апа 96 > о 


Since for rows one and two, respectively, 50 ¢ 34 and 30 ¢ 
30, all logic checks are false and the box will, therefore, 
ος ο this origin. 

Step YB. Improve the density of packing if possible. 
Since by simple inspection the box may not be moved toward 
the origin of the pallet, a detailed inspection will not be 
presented here. A detailed inspection will be given at a 


subsequent box. 
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Step 8. The second box is supported since the example 
deals with only two dimensions. 


Step 9. Update matrices A and B which are now: 


ПЕ 50 О O 60 10 ае 
А = 1 О О О 30 30 96 


2 О О О 34 o4 96 


UJ 
il 
H Hj 
t] мј 
mj "uj 


ЕЭ 
Hj 
Hj 


Step 5 Select box number 3 (20x20x96). 

me» 6. Select origin (30,0,0), B=(2,1). 

Step 7 Box will fit and improvement is not possible. 
Step 8. This box is supported. 


Step 9. Adjust À and B matrices as follows: 


1: 90 О О 60 10 96 
A = 1 О О О 30 30 96 
2 O 30 O 34 94 26 
S 30 O О 50 20 96 
F F F 
B = F Е Р 
ji ЕЁ Е 
Е Т F 


Step 5. Select box 4 (26x10). 
Step 6. Select origin (34,30,0), B=(3,1). 


Step τα. Box will fit. 
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οσο ο. mprove density. Inspect each row of matrix 
an improving direction in the x, y, and z directions. 
For illustration, the improvement in the y direction will be 
shown. It is first necessary to find the row of matrix À 
in which the following relationship does not hold (as dis- 
cussed in Section IIIE above). Note the origin is (ORX, 


ORY, ORZ) = (34,30,0) and that the box dimensions are (BOXL , р 
4 
BOXW, , BOXH, ) 

ES A 


For illustration, the inspection of row four of matrix A 


= (26x10x96). 


is shown. This corresponds to box number 3. 


AC4,2) > (34726) or A(4,3) > 30 or A(4,4) > 96 or 


A(4,5) < 34 or A(4,7) < 0 or A(4,6) < 30. 


forse results in: 30 > 60 or O > 30 or 0 > 96 от 


50 < 34 or 96 < O or 20 < 30 


Since the last term (20 < 30) is true, improvement in the 
y direction is possible. This improvement is 30-20 = 10. 
Therefore, the origin of the box is now defined to be 
(34,20,0). A similar inspection in the x direction would 
show an improvement of 4 units, making the final origin 
(30,20,0): 

Step 8. Box is supported. 


Step 9. Update matrices À and B as follows: 
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Ш 50 0 0 6ο 10 96 
A = 1 0 0 0 30 30 96 
2 0 30 O 34 54 96 
3 30 0 0 50 20 J6 


Matrix A now shows the pallet as it was loaded. Column 
one gives the identification number of the box, columns two, 
three and four give the origin of the box, and columns five, 
six and seven give the orientation of the box. 


Figure 4 shows the result of the above example. 


G. OPTIMIZATION 

In order to move from a local optimum toward the global 
optimum the following branching search can be made as 
suggested by Stoyan [4]: 

Step 1. An initial sequence AS of the boxes BOX, , i= 
(1,n), is chosen, the boxes are palletized, and the pallets 
are loaded in the container. 

πο 2 Using a uniform (1,n) pseudorandum number genera- 
tor, s random numbers are generated (s<<n) and the boxes with 
these s indexes are shuffled. 

Step 3. The new sequence is loaded and the efficiency 
of the new local optimum is measured. If the efficiency 
improved go to step number 2. 

Step 4. Replace the sequence as it existed prior to the 


shuffle and go to step 2. 
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The above algorithm is repeated until either aà predeter- 
mined amount of time is consumed or a predetermined efficiency 


is reached. 


ΕΡΤ A BETTER INITIAL SOLUTION 

Several approaches were found useful in obtaining a better 
initial solution. These included defining the input sequence 
of the boxes to be loaded according to a preconceived routine, 
selecting an optimal maximum number of turns of the input 
boxes, and sorting the pallets prior to stuffing the pallets 
into the container. These approaches are discussed at the end 


of the next chapter. 








IV. VERIFICATION 


A. RELATED ALGORITHMS 

Because the problems solved by Galata and Stoyan [4] and 
by DeSha [2] are subsets of the stuffing problem, verifica- 
tion of the current stuffing algorithm was possible by (a) 
solving the same example as presented by Galata and Stoyan 
and (b) by using DeSha's FORTRAN program to solve a problem 
which had been also solved by the stuffing algorithm. 

A solution to the minimization example of Galata and 
Stoyan was obtained by the stuffing algorithm in 0.9 seconds 
with a minimum value of z = 23.2 as opposed to the original 
solution of z = 24.5 obtained by 540 searches which took 
approximately 31 seconds each. According to reference 4, 
in their problem Galata and Stoyan loaded 95 rectangular 
solid items onto a base with the objective of minimizing 
the height, z, which was needed in order to fit all 95 
solids onto the base. In their problem, boxes need not be 
Supported. 

The FORTRAN program given by DeSha was used to solve 
the sample data discussed below: DeSha's solution showed a 
86.9% volume efficiency and an area efficiency of 92%. The 
stuffing algorithm showed an 89.1% volume efficiency and a 


95% area efficiency. 
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B. SAMPLE DATA 

mMm order to Lully test the stuffing algorithm, actual 
data were collected at Naval Supply Center, Oakland, Calif- 
ornia, Navy Exchange Retail Distribution Center in December 
1978. These data, as shown in table I, were gathered by 
individually measuring boxes which were actually loaded 
(stuffed) into an eight-foot by eight-foot by forty-foot 
shipping container. The container was loaded in six hours 
by two men who utilized one forklift truck. The men were 
asked to load the container as efficiently as possible in 
order to measure their abilities against simulated stuffing 
by a computer program. The cargo was not palletized ex- 
ternally to the container but pallets were used implicitly 
in the container. The pallets consisted of larger boxes 
placed on the floor of the container upon which smaller 
boxes were stacked. All boxes which were loaded were load- 
bearing boxes, and weight and center of gravity considera- 
tions were not addressed by the loading crew. Actual effec- 
tiveness achieved by the loading crew was 87% as compared 
to reported Naval Supply Center, Oakland averages of 90% 
for Navy Exchange cargo and 80% for general cargo (which 
consists of repair parts, equipment, and general consumable 


Supplies). 


C. SAMPLE DATA RESULTS 
The above sample data were utilized to test the stuffing 
Jeor utilizing various ordering of the input stream of 


boxes and various number of turns allowed. Table II gives 
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the results of a nine by seven factorial experiment layout. 
During this experiment, pallets were defined to be the next 
box in the ΕΓ input stream of boxes to be loaded. Each 
cell of table II gives the numbers of pallets necessary to 
stack all boxes, the IBM 360 computer run time in seconds 
necessary to load the boxes onto pallets and to stuff the 
pallets into containers, and the percent volumetric effi- 
ciency realized. Table III performs the same experiment ex- 
cept pallets were defined to be a standard pallet (40 inch 
by 44 inch) if the next box in the input. stream of boxes 
could be contained in this standard pallet. If the box could 
not be contained, the pallet was defined to be the box itself. 

Sample output of the FORTRAN program is presented in 
tables IV through VII for the test case in which standard 
pallets were not utilized and the input stream of boxes was 
ordered (highest to lowest) according to the boxes' base 
perimeter (sort control 4 as defined in the FORTRAN program) 
shown in the Computer Program section below. Table IV shows 
the boxes after they were ordered, table V shows the first 
pallet which was stacked, table VI presents a summary of 
all pallets which were stacked, and table VII shows the 
results of stuffing the pallets into the container. 

Similar results are shown in tables VIII through XI for 
the test case in which standard pallets were utilized and 
the input stream of boxes was ordered according to the boxes’ 


height. 





D. ANALYSIS OF VARIANCE 
An analysis of variance of tables II and III was per- 
formed in order to determine if a difference in efficiency 
of loading existed due to different methods of sorts and 
due to different number of turns allowed of the box. A 
level of significance of a =.10 was chosen. These analyses 
Showed that for the case in which standard pallets were not 
utilized, there was a significance between the type of sort 
utilized and there was a difference between the number of 
turns allowed. (The F statistic for these tests was, re- 
= 3.78 and F 


spectively, = 57.9.) These 


Ἵν» (8,48) 
analyses for the case in which standard pallets were utilized 
showed there was a difference in the mumber of turns allowed 


CF = 1.97) and that there was a difference in the type 


(6,48) 


o Sort utilized (F = IL) 


(8,48) 
Although the data in tables II and III indicate that 
strong interaction may exist between the number of turns 
allowed and the type of sort utilized, the lack of replica- 
tions of the experiment due to the existence of only one set 


Of data make the existence of possible interaction impossible 


to verify. 


E. RANGE TEST AND CONFIDENCE LIMITS 

In order to determine which type of sorts in tables Il 
τ 1; ава азса number of turns in table II and table III 
produced the highest efficiences, Newman-Keuls range tests 


[13] were performed. The results of these tests are shown 





55 ΧΙΙ and XIII in which an 'x' indicates a difference 
in the means and an '0' indicates no difference in the means. 
Of course, the range test was performed at the level of signi- 
meeance Of a = .10. 

Tables XII and XIII also give 90-percent confidence limits 


on the means listed therein. 


ENEEDESCUSSION OF TECHNLQUES TO IMPROVE SOLUTIONS 

Table XII shows that the greatest efficiencies when load- 
ing boxes without standard pallets occurs when boxes are 
sorted by area prior to their stacking. Stoyan's procedure 
[4], discussed in Chapter V, was then applied in an attempt 
to improve on that solution. No improvements were found in 97 
separate trials with four pairs of boxes interchanged on each 
trial. These trials tend to confirm that for the specific 
data available and for the specific method selected to load 
the boxes, the area ordering represents the best initial local 
optimum for the stuffing algorithm (when loading boxes without 
standard pallets). A similar test was conducted in the case 
where standard pallets were utilized. In this test 12 sepa- 
rate trials were conducted during which no improvement was 


Found. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The stuffing algorithm presented in this paper has 
demonstrated the possibility of achieving slightly better 
"loading" performance than is usually obtainable by exper- 
ienced loading personnel. However, its major advantages 
may be in its ability to allow reasonably accurate predic- 
tions of container requirements and in the increased speed 
with which a container may be loaded because the loading 
personnel have a plan to follow. 

This type of algorithm may be capable of allowing the 
full potential of mechanized warehouses to be realized by 
allowing the shipping department to "call forth" issues 
from the storage department when transportation assets are 
available. This would allow issue documents to accumulate 
in the mechanized warehouse's computerized data system as 
opposed to accumulating the issued material on the shipping 
dock. 

Additional data are required to verify or disprove the 
results presented in this study and to determine the exist- 
ence of the suspected interaction in the analysis of variance 
for the sorting of the boxes and the effects of the number 


of turns allowed. 
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SUMMARY LINE FOR EACH PALET THAT WAS STACKED 
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SUMMARY OF ALL PALLETS LOADED (LOADING WITHOUT STANDARD PALLETS) 
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(LOADING WITHOLT STANDARD PALLETS) 


SUMMARY GF ALL PALLETS LOADED 
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SUMMARY OF ALL PALLETS LOADED (LOADING WITHOUT STANCARC PALLETS) 
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SUMMARY OF ALL PALLETS LOADEO (LOACING WITHOUT STANDARD PALLETS) 
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